METHOD FOR LOCATING DIFFICULT ACCESS POINTS ON A MAP 

The present invention pertains to the locating of 
difficult access points, on a topological map plotted 
5 on the basis of a map of curvilinear distances. 



When dealing with a map of the zone overflown by an 
aircraft, plotted on the basis of a map of curvilinear 
distances taking account of the vertical flight profile 

10 of the aircraft, the difficult access points, which are 
those whose curvilinear distances greatly exceed the 
Euclidean distances, correspond to relief zones that 
are dangerous for the aircraft, the description 
dangerous applying to any relief zone that cannot be 

15 crossed directly by the aircraft starting from its 
current position having regard to its turning and 
climbing performance. 

The applicant has already proposed, in a French patent 

20 application filed on 9/26/2003, under no. 0311320, a 
method of estimating, on a map extracted from a terrain 
elevation database, curvilinear distances separating 
the points of the map, from a reference point taken as 
origin of the distances having regard to obstacles to 

25 be detoured around, the contours of which may change in 
the course of the time of traversal of the curvilinear 
distances as is the case for an aircraft whose current 
position corresponds to that of the point taken as 
origin of the measurements of the distances and which 

30 has to comply with a vertical flight profile with 
variations in altitude implying that one and the same 
relief that is threatening at a certain moment is no 
longer so at another or vice versa. This method 
implements a propagation-based distance transform also 

35 known by the name of chamfer mask distance transform 
since it uses a so-called "chamfer mask" array 
cataloging the approximate values of the Euclidean 
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distances separating a point of the map from its 
nearest neighbors. 



The array formed by the curvilinear distances estimated 
5 for the set of points of a map is called, for 
convenience, a map of curvilinear distances. It is not 
particularly intended to be displayed but rather to 
serve in the plotting of maps to be displayed showing 
certain specifics of the relief. 

10 

In the case of an aircraft, the map of curvilinear 
distances relates to the region overflown and has, as 
reference point taken as origin of the measurements of 
the curvilinear distances, a point near the current 

15 position of the aircraft. It serves for the plotting of 
a map, often in two dimensions, which is displayed on 
the instrument panel and shows, in false colors, a 
split of the region overflown into zones delimited as a 
function of the capacity of the aircraft to cross them 

20 and of the time that the latter would take to reach 
them when they are crossable, for example red for 
uncrossable reliefs, no route being possible, yellow 
for reliefs that are far away or close in the sense of 
the Euclidean distance but are only crossable by a 

25 diverted route and green for reliefs that are close in 
the sense of the Euclidean distance, and are crossable 
by a direct route. 

A map of the relief overflown, established on the basis 
30 of a map of curvilinear distances has the drawback of 
not giving very explicit information on the importance 
of the diversion to be accomplished when it is 
necessary to make one, thereby prompting us to 
understate, through caution, the zones represented in 
35 yellow in favor of those represented in red. 



It is possible to obtain this information on the 
importance of the diversion to be accomplished, on the 
basis of the calculation of the Euclidean distances and 
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of their comparisons with the curvilinear distances but 
account has to be taken in these comparisons of the 
presence of the obstacles to be detoured around and 
this leads to a considerable increase in the 
5 calculations required for the plotting of the map 
displayed. 

The present invention is aimed at combating this 
drawback, by depicting, on a relief map, established on 

10 the basis of a map of curvilinear distances, graphical 
information on the importance of the diversion required 
to access a point and hence, for an aircraft, on the 
dangerousness of the relief at this point, without 
however calling explicitly upon the calculation of the 

15 Euclidean distances. 

Its subject is a method of locating difficult access 
points on a topological map established on the basis of 
a map of curvilinear distances, which is noteworthy in 

20 that the map of curvilinear distances is analyzed by 
means of a chamfer mask cataloging the approximate 
values of the Euclidean distances separating a point of 
the map from its nearest neighbors, so as to extract 
therefrom, at each point of the map of curvilinear 

25 distances, the discrepancies of curvilinear distances 
separating the point considered from its nearest 
neighbors, compare these discrepancies with the 
approximate values of the Euclidean distances of the 
chamfer mask and describe the point considered as 

30 difficult of access when a difference appears. 

Advantageously, the difference noted is compared with 
several thresholds so as to devise degrees in the 
description as difficult of access. 

35 

Advantageously, the points of the map of curvilinear 
distances that are regarded as difficult of access are 
located on the topological map established on the basis 



of the map of curvilinear distances by a pattern and/or 
a particular texture. 

Advantageously, when several comparison thresholds are 
5 used to devise degrees in the description as difficult 
of access, these degrees are evidenced on the 
topological map by different patterns and/or textures. 

Advantageously, the chamfer mask used for the locating 
10 of the difficult access points is of dimension 3x3. 

Advantageously, the chamfer mask used for the locating 
of the difficult access points is of dimension 5x5. 

15 Other characteristics and advantages of the invention 
will emerge from the description below, of an exemplary 
embodiment. This description will be offered in 
conjunction with the drawing in which: 

a figure 1 represents an exemplary map of 
20 curvilinear distances covering a zone in which 

a craft is deploying and having the position of 
the craft as origin of the distance 
measurements, 

a figure 2 represents an exemplary chamfer mask 
25 usable by a propagation-based distance 

transform, 

figures 3a and 3b show the cells of the chamfer 
mask illustrated in figure 2, which are used in 
a scan pass in lexicographic order and in a 
30 scan pass in inverse lexicographic order, 

- a figure 4 illustrates the concept of direct 
trajectory for an aircraft, 

- figures 5a, 5b and 6a, 6b illustrate, as 
vertical and horizontal projections, a flight 

35 situation in which a relief constitutes an 

obstacle uncrossable by the shortest trajectory 
but crossable by a detour trajectory, 
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- a figure 7 shows the flight profile adopted for 
the map of curvilinear distances, shown in 
figure 1, 

a figure 8 shows the vertical and horizontal 
5 profiles of a relief configuration 

corresponding to a particular zone of the map 
of curvilinear distances of figure 1, 
exhibiting a partially uncrossable edge (11)/ 
a figure 9 shows an indexation used for the 
10 individual locating of the elements of the 

chamfer mask of figure 2, and 

a figure 10 is a logic chart illustrating the 
main steps of an analysis , done in a method of 
locating according to the invention, by means 
15 of a chamfer mask. 

A map of distances over a zone of deployment is made up 
of the whole set of values of the distances of the 
points placed at the nodes of a regular mesh of the 

20 zone of deployment with respect to a point of the zone, 
taken as origin of the distance measurements. As shown 
in figure 1, it may be presented in the form of an 
array of values whose boxes correspond to a splitting 
of the zone of deployment into cells centered on the 

25 nodes of the mesh. The regular mesh adopted is often 
that of the points of a terrain elevation database 
covering the zone of deployment. When a map of 
distances is used for the navigation of a craft, the 
zone point taken as origin of the distance measurements 

30 is the node of the mesh closest to the projection on 
the ground of the instantaneous position of the craft. 

Maps of distances are often produced using a 
propagation-based distance transform also known as a 
35 chamfer mask distance transform. 

Chamfer mask distance transforms appeared initially in 
image analysis to estimate distances between objects. 
Gunilla Borgefors describes examples thereof in her 
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article entitled "Distance Transformation in Digital 
Images" published in the journal: Computer Vision, 
Graphics and Image Processing, vol. 34, pp. 344-378 in 
February 1986. 

5 

The distance between two points of a surface is the 
minimum length of all the possible routes over the 
surface starting from one of the points and finishing 
at the other. In an image formed of pixels distributed 

10 according to a regular mesh of rows, columns and 
diagonals, a propagation-based distance transform 
estimates the distance of a pixel termed "goal" pixel 
with respect to a pixel termed "source" pixel by 
constructing progressively, starting from the source 

15 pixel, the shortest possible path following the mesh of 
pixels and finishing at the goal pixel, being aided by 
the distances found for the image pixels already 
analyzed and an array termed a chamfer mask cataloging 
the values of the distances between a pixel and its 

20 close neighbors. 

As shown in figure 2, a chamfer mask takes the form of 
an array with an arrangement of boxes reproducing the 
pattern of a pixel surrounded by its close neighbors. 

25 At the center of the pattern, a box assigned the value 
0 labels the pixel taken as origin of the distances 
cataloged in the array. Around this central box are 
clustered peripheral boxes filled with non-zero 
proximity distance values and mimicking the arrangement 

30 of the pixels of the neighborhood of a pixel assumed to 
occupy the central box. The proximity distance value 
appearing in a peripheral box is that of the distance 
separating a pixel occupying the position of the 
peripheral box concerned, from a pixel occupying the 

35 position of the central box. It is noted that the 
proximity distance values are distributed as concentric 
circles. A first circle of four boxes corresponding to 
the four pixels of first rank that are closest to the 
pixel of the central box that are placed either on the 



same row or on the same column are assigned a proximity 
distance value Dl . A second circle of four boxes 
corresponding to the four pixels of second rank that 
are the pixels closest to the pixel of the central box 
5 that are placed on the diagonals are assigned a 
proximity distance value D2 . A third circle of eight 
boxes corresponding to the eight pixels of third rank 
that are closest to the pixel of the central box while 
yet remaining outside the row, the column and the 
10 diagonals occupied by the pixel of the central box are 
assigned a proximity distance value D3 . 

The chamfer mask can cover a neighborhood of greater or 
lesser extent of the pixel of the central box by 

15 cataloging the values of the proximity distances of a 
greater or lesser number of concentric circles of 
pixels of the neighborhood. It may be reduced to the 
first two circles formed by the pixels of the 
neighborhood of a pixel occupying the central box as in 

20 the exemplary distance maps of figures 1 or be extended 
beyond the first three circles formed by the pixels of 
the neighborhood of the pixel of the central box. It is 
customary to stop at first three circles as for the 
chamfer mask shown in figure 2. It is only for the sake 

25 of simplification that one stops at the first two 
circles for the map of distances of figure 1. 

The values of the proximity distances Dl, D2, D3 which 
correspond to Euclidean distances are expressed in a 

30 scale whose multiplicative factor permits the use of 
integers at the cost of a certain approximation. Thus, 
G. Borgefors adopts a scale corresponding to a 
multiplicative factor of 3 or 5. In the case of a 
chamfer mask retaining the first two circles of values 

35 of proximity distance, hence of dimensions 3x3, 
G. Borgefors gives the value 3 to the first proximity 
distance Dl which corresponds to an echelon in abscissa 
or in ordinates and also to the scale multiplicative 
factor, and the value 4 to the second proximity 
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distance which corresponds to the root of the sum of 
the squares of the echelons with abscissa and with 
ordinate ^x 2 + y 2 . In the case of a chamfer mask 
retaining the first three circles, hence of dimensions 
5 5x5, she gives the value 7, which is an approximation 
Of 5a/2 to the distance Dl which corresponds to the 
scale multiplicative factor, and the value 11, which is 
an approximation of 5^5, to the distance D3 . 

The progressive construction of the shortest possible 
path going to a goal pixel, starting from a source 
pixel and following the mesh of pixels is done by 
regular scanning of the pixels of the image by means of 
the chamfer mask. 

Initially, the pixels of the image are assigned an 
infinite distance value, in fact a number high enough 
to exceed all the values of the distances that are 
measurable in the image, with the exception of the 
source pixel which is assigned a zero distance value. 
Then the initial distance values assigned to the goal 
points are updated in the course of the scan of the 
image by the chamfer mask, an update consisting in 
replacing a distance value allocated to a goal point 
with a new lesser value resulting from a distance 
estimate made on the occasion of a new application of 
the chamfer mask to the goal point considered. 

An estimation of distance by application of the chamfer 
30 mask to a goal pixel consists in cataloging all the 
paths going from this goal pixel to the source pixel 
and passing through a pixel of the neighborhood of the 
goal pixel whose distance has already been estimated in 
the course of the same scan, in searching from among 
35 the paths cataloged for the shortest path or paths and 
in adopting the length of the shortest path or paths as 
distance estimate. This is done by placing the goal 
pixel whose distance it is desired to estimate in the 
central box of the chamfer mask, while selecting the 



10 



15 



20 



25 



- 9 - 

peripheral boxes of the chamfer mask corresponding to 
pixels of the neighborhood whose distance has just been 
updated, while calculating the lengths of the shortest 
paths connecting the pixel to be updated to the source 
5 pixel while passing through one of the selected pixels 
of the neighborhood, by addition of the distance value 
assigned to the pixel of the neighborhood concerned and 
of the proximity distance value given by the chamfer 
mask, and in adopting, as distance estimate, the 
10 minimum of the path length values obtained and of the 
old distance value assigned to the pixel undergoing 
analysis . 

At the level of a pixel under analysis by the chamfer 
mask, the progressive search for the shortest possible 
paths starting from a source pixel and going to the 
various goal pixels of the image gives rise to a 
phenomenon of propagation in directions of the pixels 
which are the nearest neighbors of the pixel under 
analysis and whose distances are cataloged in the 
chamfer mask. In the case of a regular distribution of 
the pixels of the image, the directions of the nearest 
neighbors of a pixel not varying are considered as 
propagation axes of the chamfer mask distance 
transform. 

The order of scanning of the pixels of the image 
influences the reliability of the distance estimates 
and of their updates since the paths taken into account 
30 depend thereon. In fact, it is subject to a regularity 
constraint which implies that if the pixels of the 
image are labeled in lexicographic order (pixels ranked 
in row-by-row ascending order starting from the top of 
the image and progressing toward the bottom of the 
35 image, and from left to right within a row) , and if a 
pixel p has been analyzed before a pixel q then a pixel 
p+x must be analyzed before the pixel q+x. The 
lexicographic order, inverse lexicographic order 
(scanning of the pixels of the image row-by-row from 
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bottom to top and, within a row, from right to left), 
transposed lexicographic order (scanning of the pixels 
of the image column-by-column from left to right and, 
within a column, from top to bottom) , inverse 
5 transposed lexicographic order (scanning of the pixels 
by columns from right to left and, within a column, 
from bottom to top) satisfy this regularity condition 
and more generally all scans in which the rows and 
columns are scanned from right to left or from left to 
10 right. G. Borgefors advocates a double scan of the 
pixels of the image, once in lexicographic order and 
another time in inverse lexicographic order. 

Figure 3a shows, in the case of a scan pass in 
lexicographic order going from the upper left corner to 
the lower right corner of the image, the boxes of the 
chamfer mask of figure 2 that are used to catalog the 
paths going from a goal pixel placed on the central box 
(box indexed by 0) to the source pixel, passing through 
a pixel of the neighborhood whose distance has already 
formed the subject of an estimate in the course of the 
same scan. These boxes are eight in number, arranged in 
the upper left part of the chamfer mask. There are 
therefore eight paths cataloged for the search for the 
shortest whose length is taken as estimate of the 
distance . 

Figure 3b shows, in the case of a scan pass in inverse 
lexicographic order going from the lower right corner 
30 to the upper left corner of the image, the boxes of the 
chamfer mask of figure 2 that are used to catalog the 
paths going from a goal pixel placed on the central box 
(box indexed by 0) to the source pixel, passing through 
a pixel of the neighborhood whose distance has already 
35 formed the subject of an estimate in the course of the 
same scan. These boxes are complementary to those of 
figure 3a. They are also eight in number but arranged 
in the lower right part of the chamfer mask. There are 
therefore eight paths cataloged for the search for the 
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shortest whose length is taken as estimate of the 
distance . 

The propagation-based distance transform whose 
5 principle has just been recalled briefly was designed 
originally for the analysis of the positioning of 
objects in an image but it was soon applied to the 
estimation of the distances on a relief map extracted 
from a terrain elevation database with regular meshing 

10 of the terrestrial surface. Specifically, such a map is 
not furnished explicitly with a metric since it is 
plotted on the basis of the altitudes of the points of 
the mesh of the terrain elevation database of the zone 
represented. In this context, the propagation-based 

15 distance transform is applied to an image whose pixels 
are the elements of the terrain elevation database 
belonging to the map, that is to say, altitude values 
associated with the latitude, longitude geographical 
coordinates of the nodes of the mesh where they have 

20 been measured, ranked, as on the map, by increasing or 
decreasing latitude and longitude according to an array 
with two coordinate dimensions, latitude and longitude. 

For terrain navigation of mobile objects such as 
25 robots, the chamfer mask distance transform is used to 
estimate curvilinear distances taking account of zones 
which are uncrossable because of their craggy 
configurations. To do this, a forbidden-zone marker is 
associated with the elements of the terrain elevation 
30 database featuring in the map. It signals, when it is 
activated, an uncrossable or forbidden zone and 
prohibits any update other than an initialization, of 
the distance estimation made by the chamfer mask 
distance transform. 

35 

In the case of an aircraft, the configuration of the 
uncrossable zones evolves as a function of the altitude 
imposed thereon by the vertical profile of the 
trajectory adopted in its flight plan. During the 
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formulation of a map of curvilinear distances covering 
the region overflown, this is manifested as an 
evolution of the configuration of the uncrossable zones 
during the plotting of the shortest routes whose 
5 lengths serve as estimations for the curvilinear 
distances. This evolution, during the plotting, of the 
configuration of the uncrossable zones may lead to 
sizeable discrepancies between the estimations of 
curvilinear distances made for geographically close 
10 points. 

To understand this phenomenon, it is necessary to 
recall the concept of the shortest trajectory for an 
aircraft. As shown in figure 4, a shortest trajectory 
for an aircraft seeking to reach, from its current 
position 20, an aim point 21, consists, in the 
horizontal plane: 

of a rectilinear segment 22 related to the 
inertia of the aircraft, when banking into a 
turn so as to steer toward the aim point 21, 
of an arc of a cycloid 23 corresponding to the 
turning of the aircraft pushed by the crosswind 
until it reaches the azimuth of the aim point, 
and 

of a rectilinear segment 24 between the exit 
from the turn and the aim point 21. 

In the vertical plane, the shortest trajectory is 
contingent on the climb and descent capabilities of 
30 aircraft as well on the imposed altitudes. 

Certain reliefs that cannot be crossed by a shortest 
trajectory can nevertheless be crossed by a detour 
trajectory. Figures 5a, 5b and 6a, 6b give an example 
35 thereof. 

The same relief is shown in vertical cross sections, 
according to the profile of the shortest trajectory in 
figure 5a and according to the profile of a detour 
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trajectory in figure 6a, and in horizontal projections 
in figures 5b and 6a, under the guise of two strata 30, 
31 or 30', 31. Figures 5a and 5b show an aircraft in a 
current position 32 such that its shortest trajectory, 
5 located by its horizontal projection 33 and vertical 
projection 34, intercepts the relief at 35 at the 
common boundary of the strata 30, 31. Figures 6a and 6b 
show that the aircraft, in the same current position 32 
and in the same flight configuration, nevertheless has 
10 a possibility of crossing the relief illustrated by a 
first stratum 30 1 that is higher than previously 30 and 
by the same second stratum 31, by following a detour 
trajectory shown in horizontal projection 36 and in 
vertical projection 37. 

15 

A map of curvilinear distances formulated with a view 
to aiding the navigation of an aircraft takes account 
at one and the same time of the uncrossable reliefs and 
of those only crossable by detour trajectories when, in 

20 the course of the estimations of the curvilinear 
distances, the configuration of the uncrossable zones 
is made to depend on the instantaneous altitude which 
would be reached by the aircraft along the various 
routes tested assuming that it complies with an imposed 

25 vertical flight profile corresponding for example to 
that of its flight plan. Figure 1 gives a simplified 
example of such a map of curvilinear distances 
established for aiding the navigation of an aircraft 
having a vertical flight profile in accordance with 

30 that of figure 7, that is to say having a positive rate 
of climb FPA C , as is the case for an aircraft after 
takeoff. It has been formulated with the aid of the 
simplest of the distance transforms proposed by 
Gunilla Borgefors using a chamfer mask of dimension 

35 3 x 3 with two neighborhood distances 3, 4. The 
aircraft is assumed to be at the point S and to be 
moving in the sense of the arrow. The overfly zone 
covered exhibits two reliefs that are uncrossable by 
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the aircraft, one 10 completely uncrossable and the 
other 11 only crossable by detour trajectories. 

The fact that the first relief 10 is considered to be 
5 completely uncrossable amounts to admitting that the 
aircraft never reaches a sufficient altitude on the 
various routes tested for the estimations of 
curvilinear distances. Hence, its contour does not vary 
during the plotting of the various routes tested and 
10 its points keep the infinite value of curvilinear 
distance which was assigned to them on initialization. 

The second relief 11 is assumed to have the horizontal 
110 and vertical 120 contours shown in figure 8. Its 

15 vertical profile 120 approximates that of a corner, 
with a high and sheer front edge 121, for example a 
line of cliffs, facing in the direction of the current 
position S of the aircraft and leading via a descending 
line of peaks 122 to a markedly lower rear edge 123. 

20 Its high front edge 121 facing toward the current 
position S of the aircraft is crossable only on 
condition that the aircraft has gained sufficient 
altitude. This is not the case for the shortest 
trajectory which follows the propagation axes of the 

25 chamfer mask transform having as origin the current 
position S of the aircraft and going in the directions 
of the front edge 121 of this second relief 11. On the 
other hand, the aircraft will have sufficient altitude 
to cross this second relief 11, if it has taken the 

30 time to detour round it via the rear. When traversing 
the shortest routes along the second relief 11, the 
contour of this second relief 11 narrows at the rear 
until it peters out so that the chamfer mask distance 
transform ends up finding routes that are practicable 

35 for all the points belonging to the second relief 11 
which get assigned estimations of curvilinear distances 
that are lower than the initialization value. 



A map of curvilinear distances such as that shown in 
figure 1 may serve as basis for the display of a map of 
the region overflown depicting lines of equal 
curvilinear distance forming a sort of rosette around 
5 the current position of the aircraft and totally 
uncrossable terrain contours. Through the deformations 
of the rosette formed by the lines of equal curvilinear 
distance, this map also depicts terrain outlines that 
are dangerous since they are uncrossable by a shortest 
10 trajectory but these deformations are difficult to 
interpret by looking at them. 

In order to make these dangerous terrain outlines stand 
out better, although without undertaking complicated 

15 calculations, it is proposed that use be made of the 
discontinuities between curvilinear distances of 
neighboring points. The discontinuities of curvilinear 
distance between neighboring points are detected by 
scanning the points of the map of the curvilinear 

20 distances, by means of a chamfer mask cataloging the 
approximate values of the Euclidean distances 
separating a point of the map of curvilinear distances 
from its nearest neighbors. In the course of the scan, 
each point of the map of curvilinear distances is 

25 subjected to an analysis by the chamfer mask consisting 
in charting the discrepancies of curvilinear distances 
separating the point under analysis from its nearest 
neighbors, in comparing these discrepancies with the 
approximate values of the corresponding Euclidean 

30 distances of the chamfer mask and in describing the 
point under analysis as difficult of access when a 
difference is noted between Euclidean distances and 
discrepancies of curvilinear distances. 

35 The chamfer mask used for the detection of the 
discontinuities of curvilinear distances between 
neighboring points can be of any dimensions. It is 
advantageously of dimensions 3 x 3 or 5 x 5 . 



- 16 - 

Figure 9 shows the points of the neighborhood involved 
during an analysis by a chamfer mask of dimension 
3x3. These points are the four neighbors C 0 -i, C 0 i, 
C-io, C-io nearest to the point under analysis C 0 o, either 
5 in the same row, or in the same column, the four 
neighbors C-i-i, C u , C-n, Ci-i nearest to the point under 
analysis Coo of the two diagonals and the eight 
neighbors C-i- 2 , C- 2 -i, C- 2i , C-i 2 , Ci 2 , C 2 i, C 2 -i, Ci- 2 
nearest to the point under analysis C 0 o while yet 
10 remaining outside of its row, its column or its 
diagonals . 

A way of undertaking the analysis of a point by the 
chamfer mask is illustrated by the logic flowchart of 

15 figure 10. The latter consists: 

- in the course of a first step 201, in reading 
the estimated value DT(0) of the curvilinear 
distance assigned, in the map of curvilinear 
distances, to the point C 0 o under analysis, 

20 - in the course of a second step 202, in 

investigating a particular point V of the near 
neighborhood of the point C 0 o under analysis, 
preferably a point at the periphery of the 
chamfer mask, for example the point C 2 i, 

25 - in the course of a third step 203, in reading 

the value C(V) of the Euclidean distance 
separating, according to the chamfer mask, the 
point V under investigation, from the point 
under analysis Cqq, 

30 - in the course of a fourth step 204, in reading 

the estimated value DT(V) of the curvilinear 
distance assigned, in the map of curvilinear 
distances, to the point V under investigation, 
in the course of a fifth step 205, in comparing 

35 the absolute value of the discrepancy between 

the estimated values DT(0) and DT(V) of the 
curvilinear distances read in the first 201 and 
fourth 204 steps with the value of Euclidean 
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distance C(V) read in the third step 203 so as 

to note whether or not there is equality, 

in the course of a sixth step 206, in signaling 

a difficulty of access and changing the point 

Coo under analysis if the comparison of the 

fourth step 204 culminates in noting an 

inequality, 

in the course of a seventh step 207 alternative 
to the sixth step 206 should equality be noted 
at the end of the fourth step 204, in testing 
whether all the points of the near neighborhood 
of the point C 0 o undergoing analysis and 
cataloged in the chamfer mask have been 
investigated, 

in the course of an eighth step 208, in not 
detecting any discontinuity for the point 
analyzed C and in changing analyzed point Coo if 
all the points V of its near neighborhood, that 
are cataloged in the chamfer mask, have been 
investigated, 

- in the course of a ninth step 209, in changing 
investigated point V and in looping back to the 
third step 203 if all the points V of the near 
neighborhood of the point C 0 o undergoing 
analysis, that are located in the chamfer mask, 
have not been investigated. 

The test of end of investigation of all the points of 
the near neighborhood, that are cataloged by the 
chamfer mask performed in the seventh step 207, may be 
done on the maximum value of an auxiliary index for 
enumerating these points which may still be selected in 
turn, in the same order, commencing with the ones 
furthest away for which the probability of 
discontinuity is largest and ending with the ones that 
are nearest. This order of selection is for example, 
borrowing the indexation of figure 9: C-21, C-12, C12, 
C2if C2-1, C1-2, C-1-2, C-2-1, C-1-1, C-11, Cn, C1-1, C0-1, C-10, 
Coif Cio- 
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The signaling of a difficulty of access for a point of 
the map of curvilinear distances can be done by means 
of a difficulty of access pointer associated with the 
5 estimation of the curvilinear distance and used to 
modify the appearance of the points on the map 
displayed as a function of its activated or 
nonactivated state. The difficulty of access pointer 
can present several values corresponding to several 
10 values of thresholds for the discrepancies of 
estimations of curvilinear distance separating a point 
under analysis from its nearest neighbors so as to make 
it possible to display the importance of the detours 
required by differences of pattern and/or texture. 

15 

The analysis of discontinuity of curvilinear distances 
between neighboring points emphasizes the terrain edges 
that are inaccessible by a shortest trajectory such as 
the relief 11 in figure 1 which may be shown with a 

20 particular texture or pattern on the map displayed, for 
example overscoring as at 12 in figure 1. It also 
emphasizes the contours of the terrains that are 
totally inaccessible such as the relief 10 of figure 1 
but this presents less interest, these terrains being 

25 easily locatable by the initialization value of the 
estimations of the curvilinear distances of their 
points . 



